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ISGDR- bulk or surface mode ? 
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A polarization and response function for finite systems and temperatures is derived from linearizing the Vlasov 
equation. Besides the Lindhard response function in local density approximation we obtain an additional contribu- 
tion due to the surface. This formula is applied to isoscalar giant resonances which we consider as next harmonics 
to the basic dipole mode. The spurious center of mass motion is subtracted within the polarization. The results 
based on simple Fermi liquid considerations are in reasonable agreement with experimental data. 
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The existence of isoscalar giant dipole resonance (IS- 
GDR) in nuclear matter is considered as a spurious mode 
in most text books since one associates with it a center 
of mass motion. The more surprising was the experimen- 
tal justification of a giant resonance carrying the quan- 
tum numbers of a isoscalar and dipole mode^-'^^ Conse- 
quently one has to consider higher harmonics as explana- 
tion of such a mode^-''^-'^^. Usually this mode is associated 
with a squeezing mode analogous to a sound wave^-'^^^^ 

In this letter we want to discuss the influence of sur- 
face effects on the ISGDR compression mode. We will 
show that even in the frame of the Fermi liquid model 
such mode can be understood. Moreover we claim that 
the surface effects are not negligible for reproducing the 
strength function. Consequently we give first a short 
derivation of response function including surface effects. 
This will result into a new formula on the level of temper- 
ature dependent extended Thomas Fermi approximation^^ . 

The starting point is the semiclasssical Vlasov equation 



according to 

/(p,R,i) = /o(p,R) + ^/(p,R,t) 

such that the induced density variation (5n(R, t) 
I J^^f reads 
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Here we have employed the Fouriertransform of space 
and time coordinates of (|^) to solve for 6f and inverse 
transform the momentum into the form (^. Comparing 
(^) with the definition of the polarization function (||) we 
extract with one partial integration 
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where [7'^'^* is the external perturbation and [/'"'^ the self- 
consistent meanfield. Provided wc know the response to 
the external potential without sclfconsistent meanfield, 
which is described by the polarization function 11 
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The response including meanfield x^ is given by 
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X(x,x',i) =n(x,x',t) 

/ dxidx2n(x,xi,t) ^ ' 



(5n(x2,i) 



X(x2,x',i).(3) 



Therefore we concentrate first on the calculation of the 
polarization function 11 and linearize the Vlasov equation 



With (I) and 1) we have given the polarization and re- 
sponse functions for a finite system. 

In the following we are interested in the gradient expan- 
sion since we believe that the first order gradient terms 
will bear the information about surface effects. Therefore 
we change to center of mass and difference coordinates 
R = (xi -f X2)/2, r = xi — X2 and retaining only first 
order gradients we get from after Fourier transform 
of r into q 
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where in the last equality we have assumed radial mo- 
mentum dependence of the distribution function /q. We 
recognize that besides the usual Lindhard polarization 
function as the first part of (|^) we obtain a second part 
which is expressed by a gradient in space. The first 
part corresponds to the Thomas Fermi result where we 
have to use the spatial dependence in the distribution 
functions and the second part represents the extended 
Thomas Fermi approximation. So far we did not assume 
any special form of the distribution function. Therefore 
the expression ^ is as well valid for any high tempera- 
ture polarization of finite systems. 

What remains is to show that the response function (|^) 
does not contain additional gradients. This is easily con- 
firmed by two equivalent formulations of (||), 11^ = 
1 -|- Vx and xll"^ = 1 + xV, which by adding yield the 
anticommutator 
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d{Ro - R)liJin° 
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where we assumed the density dependence carried only 
by the Fermi momentum. Now it is straight forward to 
spatially average ( |l0|) with the help of (^2|) 
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n(q,cj) = / dRR'^Il{R,q,uj) 
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= n" + n™'f. 
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Consequently the surface contribution to the polarization 
function reads finally with ( pi] ) 



[n-\x]+ = 2 + [y,x]- 



(8) 



This anticommutator does not contain any gradients up 
to second order. Therefore we have [V = (5C/'"'^/(5n] 



x(R,q,w) 



n(R,q,w) 



i-y(R,q,c^)n(R,q,c^; 



+ 0{dl). (9) 



Equation (Q) and (|7|) give the response and polarization 
functions of finite systems in first order gradient approx- 
imation. 

Now we are ready to derive approximate formulae for 
spherical nuclei. In this case we can assume q||R and we 
have 



n(i?,q,w) = n"(i?,q,w)--afl 



n°(i?,q,c^) 
(10) 



where Ho is the usual Lindhard polarization with spatial 
dependent distributions (chemical potentials, density). 
We use now further approximations. In the case of giant 
resonances we are in the regime where Imll*^ ~ uj such 
that 



no = z-imno 
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(11) 



and for small q the real part vanishes for Rell'^ ^ 1 — 
c^g^/w^. Within the local density approximation we know 
that the spatial dependence is due to the density n{R) — 
nQQ{Ro—R). Since we have for zero temperature Imll'^ oc 
Pf{n) we evaluate 
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which is real. With (|T|), ^) and (|) we obtain finally 
for the structure function 



S{q,L 



[1 - y(Ren" + ff"'-f)] + (i/imn")2 
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For small q expansion we see that the pole of the structure 
function becomes renormalized similar as known from the 
Mie mode or surface plasmon mode^^^' 
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After establishing the structure function including sur- 
face contribution we specify the model for actual cal- 
culations. We choose as mean field parameterization 
a Skyrme force following Vautherin and Brink^^ which 
leads to the isoscalar potential 
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with to = -983.4 MeV fm^, ts = 13106 MeV fm^, xq = 
0.48 at nuclear saturation density uq — 0.16 fm^'^ and the 
incompressibility of K = 318 MeV. Further we employ 
the Steinwedel-Jensen model^°^ where the basic mode in- 
side a sphere of radius Rq is given by a wave vector 
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Fig. 1. The experimental structure function (T=0) ver- 
sus theoretical values. The bulk RPA result (solid 
lines) is compared with the extended Thomas Fermi 

approximation (surface corrections, dashed hues) and 10) H- Steinwedel and J. Jensen: Z. f. Naturforschung 5 



the inclusion of collisions (dot-dashed lines) . The lat- 
ter one should be of less importance due to symmetry 
of isoscalar mode. The data suggest this case and sup- 
port surface contributions. Circles: Normalized data 
from Ref.i4) 

This would correspond to the first order isovector mode-'^^^' . 
Since this mode is spurious"'^^^ we have to consider the 
next higher harmonics^'^-' which is 
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Since the polarization function with this second order 
mode contains still contributions from the spurious mode 
we have to subtract this part'^''*) 



(20) 



In Fig.l we have plotted the experimental structure func- 
tion together with different theoretical estimates accord- 
ing to ( |2^ ) and (15). The inclusion of surface correc- 
tions (dashed lines) shifts the structure function towards 
the experimental values. The inclusion of collisions (dot- 
dashed lines), which should be of minor importance for 
isoscalar dipole mode due to cancellation of backscatter- 
ing, leads to worse results. The results support also that 
the mode is of isoscalar dipole type. 
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